Titanium Nitride - a correlated metal at the threshold of a Mott transition 
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We investigate electron correlation effects in stoichiometric Titanium Nitride (TIN) using a com- 
bination of electronic structure and many-body calculations. In a first step, the Nth-order muffin tin 
orbital technique is used to obtain parameters for the low-energy Hamiltonian in the Ti-d(t% g J-band 
manifold. The Coulomb-interaction U and the Hund's rule exchange parameter J are estimated 
using a constrained Local-Density- Approximation calculation. Finally, the many-body problem is 
solved within the framework of the Variational Cluster Approach. Comparison of our calculations 
with different spectroscopy results stresses the importance of electronic correlation in this material. 
In particular, our results naturally explain a suppression of the TiN density of states at the Fermi 
level (pseudogap) in terms of the proximity to a Mott metal-insulator transition. 



I. INTRODUCTION. 

Transition metal-nitrides have been studied for sev- 
eral decades due to their appealing properties, such 
as ultra-hardness, high melting point and high Curie- 
temperature. This combination of physical and chemical 
characteristics makes them particularly suitable for coat- 
ing applications. These materials exhibit metallic con- 
ductivity and some of them even show superconductivity, 
as for example the Nb based carbonitride 1 which has a 
transition temperature of 18K. Recent low-temperature 
transport properties of thin TiN superconducting films 2 
show a disorder-driven transition from a superconductor 
to an insulating phase in which superconducting correla- 
tions persist. These experiments performed on homoge- 
neously disordered TiN films clearly demonstrate the im- 
portant role of electronic correlations. For these reasons, 
one could expect that signatures of many-body effects 
might be also present in bulk, stoichiometric TiN. 

Electronic properties of bulk transition metal ni- 
trides have been investigated using X-ray Photoelec- 
tron Spectroscopy (XPS)^, Ultraviolet Photoemission 
Spectroscopy (UPS)&, X-ray emission^, Bremsstrahlung- 
Isochromat Spectroscopy (BIS) 7 , and Electron-Energy- 
Loss Spectroscopy (EELS)&2. A number of optical re- 
flectivity measurements have also been carried ou t 10 ' 11 . 
On the theoretical side, a large number of band struc- 
ture calculations using density functional theory (DFT) 

- mostly within the local-density-approximation (LDA) 

- is present in the literature 3 - ' 12 ! 13 ) 14 . LDA results de- 
scribe bonding in terms of (i) a metallic contribution 
giving a finite density of states at the Fermi level, (ii) 
an ionic contribution caused by the charge transfer from 
the metal to the non-metal atom, and (iii) a covalent con- 
tribution due to the interaction between the non-metal 
p and the metal d valence states in addition to the (iv) 
metal-metal interactions. In the case of TiN, it is believed 
that bonding is mostly covalent in origin. However, while 
the LDA results for the occupied part of the density of 
states (DOS) show a good agreement with spectroscopic 



data at high binding energies, they fail to describe cor- 
rectly the energy range in the vicinity of the Fermi energy 
(see, e. g., Fig. [2]) where predominantly Ti-3d states are 
presen t 15 d 6 i 17 . In particular, a suppression of the XPS 
spectrum over a range of energy down to rj leV below the 
Fermi energy, could be explained so far only by means of 
an artificial broadening, whose parameters are optimized 
to fit the experimental data. In this paper, we present an 
alternative view and argue that this suppression can be 
naturally explained in terms of correlation effects, and, 
specifically, by the proximity to a Mott metal-insulator 
transition. 

Recent developments in the methodology of electronic- 
structure calculations allow to go beyond the LDA and 
include electronic correlation effects, which are partic- 
ularly important for electrons in d or / orbitals. The 
simplest mean-field extension of LDA is the LDA+U ap- 
proach 1 -. However, as we show below, even such ex- 
tension is not sufficient to improve the agreement with 
experiments in TiN. For this reason, it is necessary to 
go beyond the mean-field approach and to appropriately 
deal with dynamical correlation effects. In this paper, 
we carry out this task by means of a cluster-perturbative 
method, the Variational Cluster Approach (VCA) 19 i 20 , 
which we combine with LDA in order to obtain the ap- 
propriate model parameters. 

This paper is organized as follows: the results of the 
electronic-structure calculations within the framework of 
the linear muffin-tin orbital method (LMTO)2i at the 
LDA level are presented in section fll] Comparison of 
LDA and LDA+U results with experiments are also dis- 
cussed in that section. In Sec. Ill Al we describe the ab- 
initio construction of the many-body model Hamiltonian. 
Specifically, the uncorrelated part of the Hamiltonian for 
excitations in the vicinity of the Fermi level is obtained 
from the so-called downfolding techniqu e 22,23 within the 
Nth-order muffin tin orbital (NMTO) method. The in- 
teraction part is then estimated by the constrained-LDA 
method. In Sec. Ill Bt we give a short summary of the 
VCA approach. We present and discuss our results in 
Sec. Mil In particular, in Sec. MI Al we evaluate the den- 



2 



sity of states within VCA and compare it with experi- 
ments, and discuss the results in the framework of previ- 
ous calculations. In Sec. IIIIBI we discuss fc-dependent 
spectral properties, namely the spectral function and 
the self-energy, and evaluate the effective electron mass 
which we also compare to experiment. Finally, we sum- 
marize our results in Sec. IIVI 



II. 



ELECTRONIC STRUCTURE 
CALCULATIONS FOR TIN 



TiN crystallizes in the rock-salt (Bl) structure where 
Ti and N atoms are sitting on interpenetrating face- 
centered cubic lattices originating at (0,0,0), and at 
(5,5,5) in units of the lattice parameter a — 7.65ao 
(do = Bohr radius), respectively. Each Ti (N) atom has 
six N (Ti) nearest neighbors in an octahedral geome- 
try. Fig. [T] shows the conventional unit cell containing 
Ti (large spheres, red) and N (small sphere, blue) atoms. 
The electronic configuration of the outer shell of Ti is 




FIG. 1: (color online) Octahedral environment formed by the 
Ti-atoms (large, red spheres) around the N atom (small, blue 
sphere). The unit cell is represented by the black lines; the 
subset of Ti-atoms connected by dashed lines shows the clus- 
ter geometry that was used as reference system for the VCA- 
calculation (see text). 

4s 2 3d 2 . Therefore, following from the formal oxidation 
Ti 3+ , there is a single Ti electron in the valence band. 
Due to the octahedral symmetry, the Ti-ci-orbitals are 
split into the three ti g orbitals at lower energy, and two 
e g orbitals at higher energy. 

The LDA band structure of TiN was computed with 
the LDA-LMTO (ASA) code 2 ^. which uses the basis of 



linearized muffin-tin orbitals in the atomic sphere ap- 
proximation. Numerically, two empty spheres per unit 
cell have to be introduced at the positions (1/4, 1/4, 1/4) 
and (3/4, 3/4, 3/4), due to the atomic sphere approxima- 
tion (ASA) used in the calculation. Results are shown 
in Fig. [5J In the density of states, the N-s orbital (not 
shown) form the lowest valence band widely separated 
from the other valence bands. At higher energies, one 
finds the set of bands formed by N-p orbitals situated 
between -10 and -4eV. Finally, the last five bands mainly 
consist of Ti-d orbitals. The DOS around the Fermi level 
is dominated by Ti-t2 g bands, while e g bands remain 
above the Fermi level and are empty. t2 g and e g bands 
overlap in the unoccupied part of the spectrum. On the 
other hand, the hybridization of Ti-d orbitals with N- 
p orbitals below the Fermi level is at the origin of the 
covalent bonding. Concerning the p — d hybridization, 
the energetically higher d-bands form the anti-bonding 
states while the lower p-bands form the bonding states. 
The position of the Fermi level (Ep) is determined by 
the number of valence electrons per unit cell being equal 
to 9. 

In order to take into account correlation effects on a 
mean field-level, we carried out an LDA+U calculation. 
Here, we used values of U = lOeV and J = 1.3eV, as 
obtained from the constrained LDA calculations^^. As 
one can see there are no significant differences between 
the LDA and the LDA+U results despite the large value 
of U. However, a comparison of these calculations with 
valence band XPS spectra, which provide a measure of 
the total density of occupied states as a function of en- 
ergy, shows that neither LDA nor LDA+U results are 
appropriate to describe the experimental spectra in the 
vicinity of the Fermi energy unless one introduces ad-hoc 
broadening terms. Specifically, while the position of the 
N-p bands obtained by LDA and LDA+U is in reasonable 
agreement with the XPS measurements (see Fig [5]), both 
methods fail to reproduce the behavior of the DOS within 
a range of ~ 2 eV below the Fermi energy. In particu- 
lar, the XPS spectra show a local maximum at energies 
of ps — feV, followed by a "pseudogap" -like suppression 
at the Fermi energy. The LDA+U results do not change 
qualitatively when decreasing U down to U « 8, which, 
as argued below, is more appropriated for this material. 

These results suggest that static inclusion of correla- 
tions is not sufficient to explain the DOS near the Fermi 
energy. For this reason, we have taken into account dy- 
namical correlation effects beyond the LDA results by 
means of the VCA, a method appropriate to treat cor- 
related systems. This approach builds up on the exact 
diagonalization of an Hamiltonian on a finite cluster com- 
bined with an appropriate extension to the infinite-lattice 
limit. However, in order to limit the size of the Hilbcrt 
space, it is necessary to use an effective Hamiltonian de- 
scribing a small number of correlated effective orbitals 
per lattice site in the close vicinity of the Fermi level. 
From Fig. [21 we conclude that the minimal model has 
to be restricted to Ti (t2 g ) bands, while higher Ti (e 9 ) 
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FIG. 2: (color online) Orbital-resolved density of states for 
Ti obtained by LDA and LDA+U. The upper (lower) panel 
shows the t2 g (e g ) contribution. The upper panel also displays 
the N-p density of states and the experimental XPS spectrum^ 
(Experimental data reproduced with kind permission of the 
authors). 



ber of energy points on which the LDA bands are re- 
produced exactly. To optimize the overall agreement 
with the energy bands, we chose the following expansion 
points E V -E F = 1.512eV, 0.016eV, and -0.664eV. Re- 
sults are quite stable for choices of the E v around these 
values. In Fig. [3] we show the eigenvalues of the effec- 
tive Hamiltonian along some high-symmetry directions 
in comparison with the bands obtained from the full or- 
bital basis. From this figure it is clear that the Ti-d(t2 g ) 
manifold is well described by the non-interacting part of 
the effective Hamiltonian. 
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bands are unoccupied and can be neglected. 



A. Ab-initio construction of the model 
Hamiltonian 

In order to construct the effective low-energy Hamilto- 
nian to use in our VCA calculation, we employed the Nth 
order muffin-tin-orbitals scheme within the same LMTO- 
ASA basis set. The NMTO method^^ can be used to 
generate truly minimal basis sets with a massive down- 
folding technique. This reduced basis set reproduces the 
bands obtained with the full basis set to great accuracy 
and can thus be used as non-interacting part of the many- 
body Hamiltonian. According to the discussion in Sec.HIl 
the minimal basis set is obtained by downfolding to the 
Ti-d(t2g) manifold. The truly minimal set of symmetri- 
cally orthonormalized NMTOs is a set of Wannier func- 
tions. In the construction of this set, the active chan- 
nels are forced to be maximally localized onto the eigen- 
channel Rim (R = atomic position, /, m — angular mo- 
mentum quantum numbers), which makes the basis set 
strongly localized and suitable for a real-space Hamilto- 
nian. 

In this way, the non-interacting part of the effective 
Hamiltonian is confined to a reduced set of orbitals and 
to a reduced energy window. The NMTO downfolded 
bands are obtained by expanding around a small num- 



FIG. 3: (color online) Band-structure of TiN calculated for 
a path in the Brillouin zone (BZ) going from the L-point 
(0.5,0.5,0.5) trough the F (0,0,0), X (0,1,0), W (0.5,1,0), 
L (0.5,0.5,0.5), K (0,0.75,0.75) points and ending at the F 
point (0,0,0). The bands are obtained from the LDA calcu- 
lation (dashed/green line) for the complete orbital basis set. 
The NMTO bands (full/red line) are obtained after downfold- 
ing to the Ti-3d(t2 9 ) orbitals. 



Fourier-transformation of the orthonormalized NMTO 
Hamiltonian, ff LDA (k), yields on-site energies and hop- 
ping integrals (R = (x, y, z)), 

(Xr/ot' \H LDA - £f| Xkm) = *<m R ' W 

in a Wannier representation, where the NMTO Wan- 
nier functions XRm) are orthonormal. In the restricted 
NMTO basis, m is labeled by the three t^g orbitals 
m = xy,yz, zx (we use this order for the definition of 
the matrix elements below). In this basis, the on-site 
matrix elements t$5?/° m are diagonal and independent of 
m. The precise value of these terms, as well as the so- 
called double-counting correction, are not important, as 
they can be absorbed into the chemical potential. The di- 
rectional hopping matrix elements up to second nearest- 
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neighbor are given by 
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in units of eV. Only one representative hopping integral 
is shown for each class. Other hopping terms can be 
derived from proper unitary transformation using crystal 
symmetry. See, e.g., Ref. 12a for details. 

Further hoppings are taken into account up to a range 
of r = 1.1a. Neglected hoppings are at least by a factor 
40 smaller than the largest nearest-neighbor hopping. 

In order to make sure that neglecting e g orbitals is 
safe, we have also carried out a LDA+VCA calculation 
(for smaller clusters) using an Hamiltonian downfolded 
to all 5 Ti-d bands. We have verified that the occupa- 
tion of the e g bands is less than 10~ 3 per Ti-atom, so 
that the corresponding Wannier functions can be safely 
neglected in our calculation. Notice, however, that, due 
to hybridization, t2 g Wannier functions also have a cer- 
tain amount of e g character. 

The non-interacting part of the effective Hamiltonian 
for TiN, thus, has the form 



u _ \ " .R'-R t 

■"0 - 2^ t m>,m C R'm'a C ^r, 

R/,R,,{m',m},fT 



(2) 



To take into account correlation effects, we add the usual 
interaction term 



ml 



(3) 



R,,m 



n 'R.ma n 'R,rn' a' 

R,,m<m / ,o",o' / 



' /c Rm'T C RmJ, C R™'l C RmT + n - C - 

Jc Rm'T c R™'l CR ™l CR ™T + h - c - 



R,m<m' 

In Eq. [31 CR, mCT (el, ) are the usual fermionic annihila- 
tion (creation) operators acting on an electron with spin 
a at site R in the orbital m and riR„ lCT = c^ mcr CR mcr . 
U denotes the Coulomb-interaction for two electrons in 
the same orbital with anti-parallel spin. If they are lo- 
cated on two different orbitals, the interaction is reduced 
to U' — U — 2 J. J is the Hund's rule exchange constant 
and h.c. denotes the hermitian conjugate. 

We have estimated the value of the Coulomb- 
interaction parameter by means of constrained LDA t 24 i 25 
whereby occupancies on all d orbitals have been fixed. 
Since the e g effective orbitals are essentially empty, as dis- 
cussed above, they do not contribute to screening and can 
be safely neglected. The constrained calculation yields 



a value of U ~ lOeV and J rs 1.3eV. However, as it 
was shown by Aryasetian et al22, constrained LDA gives 
larger values for U compared to other methods based 
on the evaluation of the screened Coulomb interaction 
within the Random Phase Approximation 28 . Therefore, 
we have also investigated smaller U values. We will show 
below that results change drastically around U ~ 9 eV, 
where a Mott-insulator transition takes place. 

As correlations on a mean-field level are already in- 
cluded in LDA, one should in principle subtract the 
long discussed double-countin g 29 i 30 correction. However, 
since ti g orbitals are degenerate, this correction is simply 
a constant that can be absorbed in the chemical poten- 
tial. 



B. Variational Cluster Approach 



To solve the many-body Hamiltonian ([3J)+ (J3J) we 
employ the Variational Cluster Approac h 19 ' 20 . This 
method is an extension of Cluster Perturbation Theory 
(CPT) 31 i 32 ' 33 , in which the original lattice is divided into 
a set of disconnected clusters and the inter-cluster hop- 
ping terms are treated perturbatively. VCA addition- 
ally includes "virtual" single-particle terms to the cluster 
Hamiltonian, yielding a so-called reference system, and 
then subtracts these terms perturbatively. The "opti- 
mal" value for these variational parameters is determined 
in the framework of the Self-energy Functional Approach 
(SFA)^! 3 ^, by requiring that the SFA grand-canonical 
potential f2 is stationary within this set of variational 
parameters. Since TiN is paramagnetic, we only include 
the chemical potential of the cluster as a variational pa- 
rameter. The latter is necessary in order to obtain a 
thermodynamically consistent particle densit y 36 ' 37 . In 
this paper, we use a new method, described in Ref. l38l . 
to carry out the sum over Matsubara frequencies required 
in the evaluation of f2, whereby an integral over a con- 
tour lying a distance A from the real axis is carried out. 
The crucial point is that the contour integral is exact for 
any (even large) A. 

As a reference system we adopt the minimal cluster 
that contains the full lattice symmetry. This consists 
of all 6 Ti sites lying on the corners of an octahedron, 
as shown in Fig. [TJ Larger clusters are at present not 
feasible within our variational procedure. As in the case 
of cluster-DMFT (which is currently not feasible for a 
6-site cluster with three orbitals each), the appropriate 
periodisation is a crucial issu o 39 ' 40 . Here, we choose to 
periodise the Green's function, as it is well known that 
the self-energy periodisation gives unphysical results in 
the vicinity of an insulating phase. 
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III. RESULTS 
A. Density of states 

As discussed in SecHH the low-energy XPS 3 spectrum 
is characterized by a peak at about —1 eV followed by a 
"pseudogap" , i. e. a suppression of states at the Fermi 
level. This latter fact is consistent with the K-ELNES 
spectrum^, which provides information about the DOS 
above the Fermi energy (cf. Fig. [5|. This suppression is 
not reproduced by LDA electronic structure calculations, 
suggesting that strong electronic correlations may play an 
important role for this material. In a previous work^, it 
was suggested, in order to improve the agreement with 
the measured spectra, to convolute the computed DOS 
with a combination of a Lorentzian and a Gaussian curve, 
taking into account the broadening due to lifetime and 
experimental resolution effects. This treatment indeed 
improves on the overall shape producing a peak at an 
energy of about —1.5 eV, although one should point out 
that the fitted broadening parameters are much too large 
(around 0.6 and 0.8eV, respectively). 

In the present paper, we argue that the pseudogap ob- 
served in the DOS of stoichiometric TiN is intrinsic to 
this material and signals the proximity to a Mott metal- 
insulator transition. To show this we start by calculating 
the DOS obtained from the Hamiltonian (J3J) with the 
value of U = 10 eV obtained from constrained LDA. The 
results, displayed in Fig. |H predict for this value of U 
a Mott-insulating state with a gap of about leV. How- 
ever, this is in contrast to the experimental situation, 
since electrical conductivity in TiN shows a metallic be- 
havior, although with a relatively low residual conduc- 
tivity^. Since results obtained from constrained LDA 
are only approximate and tend to overestimate U due to 
the delocalized nature of Wannier orbitals, we have also 
presented results for smaller values of U down to 8eV 
(Fig. 01 . As one can see, no significant changes can be 
detected for states more than 3eV away from the Fermi 
energy. Here, only static correlations, which get absorbed 
into the chemical potential, play a role. On the other 
hand, the situation changes rapidly around the Fermi 
energy. In particular, our results show that the Mott 
gap starts closing at U w 9eV, and at smaller U down to 
U k 8eV it acquires the shape of a pseudogap. As a mat- 
ter of facts, the curve for U = 8.5, shown in Fig. [5] agrees 
quite well with experimental measurements. Notice that, 
in order to avoid introducing ad-hoc parameters, and to 
show fine-detailed features of the spectrum, the calcu- 
lated curve has not been additionally broadened to meet 
experimental resolution. Taking this into account, we 
see that our results reproduce quite well the experimen- 
tal features both above (K-ELNES) and below (XPS) 
the Fermi energy, and in particular the "pseudogap" be- 
havior between —1.0 and l.OeV. We stress that we are 
using an effective low-energy model which is expected to 
correctly reproduce correlation effects close to the Fermi 
energy, but is not expected to reproduce features beyond 



the range of the NMTO bands shown in Fig. [3J (full/red 
curves). To reproduce the spectrum at higher binding 
energies, LDA and LDA+U are more appropriate. In 
this sense, our results complement these techniques in 
the low-energy region. 




E-E F (eV) 

FIG. 4: (color online) TiN density of states calculated with 
LDA+VCA for different values of the Coulomb-interaction 
parameter U and for J = 1.3. 

Another remarkable feature of the computed 
LDA+VCA density of states spectra is the pres- 
ence of a set of peaks in the energy range of —5 to 
— 3eV. These states only appear when a sufficiently large 
cluster is taken as a reference system. Therefore, these 
states can be seen as non-local many-body incoherent 
features. It is tempting to associate these states with the 
corresponding "kink" in the XPS spectrum^ at energies 
around — 4eV. However, this kink is also affected by the 
presence of N-2p states whose bands start at w —4eV 
(see Fig. [2]). Therefore, our result suggests a strong 
hybridization between these states and the correlated 
many-body structures at « — 4eV. 

The nitrogen K-ELNES spectrum is related to the un- 
occupied partial density of states with p-symmetry at the 
nitrogen site. The peaks situated at energies around 2 
and 4eV can be attributed to the unoccupied N-2p states 
hybridized with Ti 3d states. Notice that the downfolded 
operators used in our effective Hamiltonian © + ([21) de- 
scribe anyway an effective particle ultimately producing 
the bands of Fig. [3J i.e. the corresponding particle also 
contain a "mixture" of other orbitals such as N-p, in or- 
der to correctly reproduce the hybridization. It is re- 
markable to note that our calculation with U = 8.5eV 
captures most of the features of the experimental spec- 
tra at both peaks, in particular the separation of the two 
peaks is in very good agreement with the experimental 
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FIG. 5: (color online) TiN density of states calculated within 
LDA+VCA for U = 8.5eV and J = 1.3 compared with mea- 
sured XPS 3 and K-ELNES 17 spectra (Experimental data re- 
produced with kind permission from the authors). 



value 2.3eVii. 

It is important to mention that combining the exper- 
imental data from the occupied and unoccupied parts 
of the spectra with electric and magnetic properties of 
TiN 4 ^, a peculiar metallic behavior emerges. Magnetic 
susceptibility measurements show that TiN is a Pauli- 
paramagnet, and electrical conductivity demonstrates a 
metallic behavior with a relatively large residual resis- 
tivity. In combination with the measured XPS spectra, 
one can conclude that at the Fermi level a pseudogap in 
the density of states is formed, signaling the vicinity of a 
Mott-insulating phase. 

The present LDA+VCA calculation includes correla- 
tion effects exactly on a length scale of an octahedral 
cluster consisting of 6 sites shown in Fig. [TJ (connected 
by dashed lines) . Notable results of our present calcu- 
lations are the correct description of the — leV peak and 
the double peak in the occupied part of the spectra. In 
addition, at lower energies, many-body non-local inco- 
herent features are formed which can explain the kink 
in the XPS spectra situated at — 4eV. Therefore, it is 
clear that the present LDA+VCA results show notable 
improvements with respect to previous DFT-results in 
the low-energy region, and explain the peculiar metallic 
behavior of TiN. 

A remark should be made concerning vacancy effects. 
It is well known that transition metal carbides and ni- 
trides usually contain vacancies in the metalloid (N) sub- 
lattice. The presence of vacancies profoundly influences 
the physical properties of this family of compounds. It 
is known from electronic-structure calculations that the 



presence of vacancies reduces the partial nitrogen s and 
p density of states, and produces additional peaks close 
to the Fermi level. These vacancy peaks show up in the 
LDA-DOS of non-stoichiometric TiN at about — 2eV and 
in the vicinity of the Fermi leve l 15 i 16 ' 43 i 44 . The peak at 
— 2eV is associated with cr-bonding between Ti atoms 
through the N vacancy, while the peak at Ep is re- 
lated to the cr-bonding between the nearest-neighbor Ti 
atom a 15 i 43 . A quantitative comparison of the measured 
spectra and the broadened density of states obtained for 
non-stoichiometric materials shows that the calculated 
density of states peaks are too narrow and shifted to lower 
binding energies with respect to experimental results. 
The shift and the narrowing of the theoretical peaks was 
interpreted 4 ^ as due to limitations of the local density 
approximation in describing electronic interactions. 



B. Spectral properties 

In order to gain insight into k-dependent features of 
the local DOS, we have computed the k-resolved spectral 
function A(k, u>), which is plotted in Fig. [5] From this 
figure one can clearly discern two metallic bands crossing 
the Fermi energy near to the W point, and one between K 
and r, however, with a small spectral weight, consistent 
with the pseudogap picture discussed in Sec. IIII Al 
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FIG. 6: Spectral-function of TiN for U = 8.5 and J = 1.3 
shown as density plot for the same path in the BZ as in Fig. [3] 
The bands obtained from downfolding (dashed, green lines) 
are shown as reference. 

Moreover, in the energy range between —3 to — 5eV 
the spectral function in Fig. [6] shows dispersion-less fea- 
tures which are responsible for the non-local correlation 
peaks discussed in the previous section. As discussed 
above, these features cannot be captured within a single- 
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site LDA+Dynamical-Mean-Field (DMFT) approach, so 
that a cluster approach, as the one presented here, is 
required. Also the cluster geometry itself is important: 
we chose an octahedral reference cluster (Fig. [1} in order 
to conserve the lattice symmetry. This is needed, as we 
have verified that in smaller clusters, without the lattice 
symmetry, this feature is not present. 

As discussed above, our calculations show that dynam- 
ical correlations are important for TiN. On the other 
hand, it turns out that also non-local correlations are 
crucial. An important consequence of this fact is that we 
do not expect the properties of TiN discussed here to be 
correctly reproduced by a single-site DMFT calculation. 

To show this, we plot in Fig. [7] the self-energy on the 
same path around the BZ as for the spectral function 
(Fig. [6]). As one can see, there is a strong k dependence, 
especially in the region around the Fermi energy. This is 
an indication for the non-locality of the self-energy. This 
non-locality strongly affects the metal-insulator transi- 
tion as well, since the most dispersive part of £ is pre- 
cisely around the Fermi energy. Furthermore, there are 
contributions to the self-energy at larger energies (not 
shown in Fig. [7]), which, however, are essentially flat as 
a function of k, and, thus, localized in real space. 




E-E F (eV) 

FIG. 7: Imaginary part of the self-energy calculated for the 
same parameters and along the same BZ path as in Fig. [6] 
and presented in a three-dimensional plot. 

The k-resolved spectral function allows us to determine 
the LDA+VCA quasiparticle bands, which are shown as 
density plot in Fig. [5] and are compared with the bands 
obtained from LDA (dashed). One can observe a reduc- 
tion of the quasiparticle bandwidth (BW) from approx- 
imately 7.5 eV in LDA to 4.5 eV in LDA+VCA. Ac- 
cordingly, the band narrowing due to correlation effects 



can be described in terms of a "high-energy" (HE) mass- 
renormalization factor ttihe j ttihe.lda ~ 1-7. Experi- 
mental measurements provide the total ratio run e /m e 
(m e is the free electron mass), which contains also the 
mass renormalization factor TnnE,LDA/ m e coming from 
the band structure. The latter can be roughly esti- 
mated by equating h 2 k 2 F / {2ijihe,lda) = E* F , whereby 
E* F = AA2eV is the Ti-3d band "depth" , i. e. the Fermi 
energy measured from the bottom of the occupied Ti-3d 
band, and hp is related to the corresponding occupation 
n = 4/a 3 by the usual relation kp = 3ir 2 n. This gives 
m HE,LDA/ m e ~ 1.3, producing an overall high-energy 
mass renormalization m#£/m e = 2.2. This has to be 
distinguished from the low-energy effective mass given by 
the change in slope at the Fermi energy, i. e. the Fermi 
velocity. Nevertheless, available experimental measure- 
ments discussed below provide an estimate for muE and 
not for the low-energy contribution. 

Experimentally, an estimate of the electron effective 
mass in TiN has been obtained by UPS and EELS exper- 
iments using high-energy synchrotron radiation^. These 
methods also allow to determine the band density n of 
the ci-electrons and the 3d band depth (corresponding to 
our E F above). To estimate the (i-electron density several 
characteristics of the UPS spectrum such as peak area, 
photoemission cross section and inelastic mean free path 
are used. The estimation of the band depth is obtained 
from the position of the minimum in intensity between 
the valence and conduction bands^. Considering the 
rough estimate we have used, the computed high-energy 
mass renormalization factor m#£/m e = 2.2, turns out 
to be in reasonable agreement with the value 2.7 ± 0.3, 
obtained in Ref. [46|, table 3, for the stoichiometric TiN 
sample d. 



IV. SUMMARY 

In this paper we have analyzed the physical properties 
of TiN in a combined electronic-structure and many-body 
approach. The NMTO downfolding technique was used 
to calculate the LDA low-energy effective Hamiltonian 
for the Ti — dfcg) bands. The effect of non-local corre- 
lations was treated within a variational cluster perturba- 
tion approach, which deals exactly with correlations on 
the cluster scale. We present results for the local density 
of states, the self-energy, the spectral function and for the 
effective mass. The results are analyzed in comparison 
with experiments on valence band XPS and K-ELNES 
spectra and, in the case of the effective mass, compared 
to UPS and EELS experiments. 

Our results suggest that TiN is a peculiar metal with 
a pseudogap at the Fermi level indicating the proximity 
to a metal-insulator transition. In our calculations the 
pseudogap regime is best described for a value of U = 
8.5eV for the Coulomb-interaction. 

It is important to mention that this result could be 
achieved only by an appropriate treatment of dynam- 
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ical and non-local correlations. In particular, neither 
LDA nor LDA+U calculations provide even a qualita- 
tive description of the DOS suppression at the Fermi en- 
ergy. Due to the non- locality of the self-energy, we expect 
that even single-site DMFT may not be appropriate for 
a proper description of this system at low energies. 

Our calculations provide a good qualitative and semi- 
quantitative agreement with experiments, when compar- 
ing the DOS with XPS and K-ELNES spectra close to the 
Fermi energy. In particular, our results show a coherent 
feature situated at — leV which can be clearly identified 
in the XPS spectra. We stress that none of previous 
theoretical investigations are able to capture this partic- 
ular low energy feature. Furthermore, at low energies our 
results provide a qualitative interpretation of the — 4eV 
kink in the XPS spectra as an hybridization of incoherent 
non-local many-body features with the N-p states. Also, 



in the unoccupied part above the Fermi level, our results 
reproduce the double peak structure in good agreement 
with the K-ELNES spectra. Finally, we calculated the 
effective mass from high- and low-energy spectral fea- 
tures and found reasonable agreement with experimental 
values. 
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